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ABSTRACT 

This  paper  is  to  prove  that,  the  fractional  partial  differential  equation  forms  a  finite  domain  with  the  numerical 
solution  by  using  different  fractional  derivatives.  The  two  fractions  are  used  to  prove  the  fractional  diffusion  equation  and 
the  fractional  dispersion  equation.  The  Fractional  Differential  Equation  is  formed  from  the  standard  diffusion  equation 
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replacing  the  order,  as  a  second  space  derivative,  with  the  fractional  derivative 
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"  '  '.  The  analytical  solutions  of,  both  the  fractional  diffusion  equation  and  dispersion  equation  were  derived. 
From  these  three  numerical  methods,  the  Ll/L2-approximation  method,  the  standard  method,  and  the  matrix  transform 
method;  only  the  third  method  is  used  to  deal  with  the  fractional  derivative.  The  two  fractional  methods  transmute  into  a 
system  of  ordinary  differential  equations  that  solves  by  the  graphical  method.  It  concludes  by  the  numerical  results,  that 
have  demonstrated  the  three  numerical  method’s  effectiveness  and  convergence. 
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1.  INTRODUCTION 

Fractional  derivative  is  not  a  new  concept,  but  it  is  an  old  and  familiar  integer-order.  The  fractional 
derivatives  were  applied  in  many  fields,  like  hydrology  physics,  biology,  chemistry,  finance,  biochemistry,  etc.  This 
new  fractional- order  model  is,  more  adequate  than  the  initially  used  integer-order  models,  the  reason  is  that,  the 
fractional- order  derivatives  and  integrals  need  the  description  of  the  memory  and  hereditary  properties  of  different 
substances.  This  is  the  important  advantage  of  the  fractional- order  models,  when  we  compared  with  the  neglected 
integer- order  models. 

In  Physics,  the  classical  Brownian  motion  model,  spreads  with  a  particle  in  an  inconsistent  way,  when  the 
fractional  space  derivatives  are  merged  to,  model  anomalous  diffusion  or  dispersion.  Similarly,  a  left  Riemann- 
Liouville  derivative  and  right  Riemann-Liouville  derivative,  were  used  to  allow  the  modeling  of  flow  regime  impacts, 
from  either  side  of  the  domain.  To  model  the  transport  of  passive  tracers,  carried  by  fluid  flow  in  a  porous  medium,  for 
solute  transport  in  a  subsurface  material  dispersion  equation,  were  used  in  groundwater  hydrology.  Green  functions  are 
difficult  to  evaluate  but,  the  fractional- order  partial  differential  equations,  were  analytically  derived  by  it.  It  proves 
that,  its  great  potential  is  helped  to  derive,  the  numerical  treatment  and  it  supports  to  analyze  the  fractional- order 
differential  equations. 

Probably,  we  cannot  use  the  numerical  method,  to  solve  the  fractional  partial  differential  equation,  with  non- 
symmetric  fractional  derivatives.  The  practical  numerical  method  was  used,  to  solve  the  one- dimensional  space,  with 
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variable  coefficients,  on  a  finite  domain.  The  backward  differentiation  formula  was  also  used,  to  solve  the  space  fractional 
equation,  into  a  system  of  ordinary  differential  equations.  The  domain  decomposition  method  and  variation  iteration 
method,  are  the  two  reliable  developed  algorithms,  which  are  used  to  construct  numerical  solutions,  that  will  be  in  the  form 
of  a  rabidly  convergent  series,  with  computable  components,  so  that,  the  theoretical  analysis  cannot  be  proved  by  that. 

2.  THE  FRACTIONAL  DIFFUSION  AND  DISPERSIOIN  EQUATION 

In  this  paper,  the  functionall  partial  differential  equation  forms  a  finite  domain,  with  the  numerical  solution,  by 
using  different  fractional  derivatives. 

cfaU,  i)  &  d1* 

~  K*  t} + £) ' 

0  <  t  <T}  0  <  x  <  L. 


It  is  given  by  the  initial  boundary  as 

id.Q,  t\  =  u{L,  i)  =  0, 
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The  solute  concentration  is  ‘u’ 


js  dispersion  coefficient 

KqKq 

p  p  Is  the  average  fluid  velocity 

The  Fractional  derivatives  of  order 
defined  as 
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Different  techniques  were  used  for  different  fractional  derivatives,  without  confusing  the  remark  that  is  essential 
to  adopt  the  accuracy,  for  the  corresponding  definition  of  the  fractional  operator,  for  the  standard  accuracy  and 
discreteness,  of  the  fractional  diffusion  equation.  The  unstable  finite  difference  scheme,  results  in  finite  difference  method, 
as  an  explicit  or  an  implicit  one. 
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Fractional  Laplacian  operator  and  derivation  are  equivalent;  this  was  widely  assumed  by  the  fractional  derivative 
method.  Firstly,  the  Laplacian  was  used  by  the  Fourier  transform,  on  an  infinite  domain  with  a  natural  extension,  to  include 
finite  domains.  The  second  is  the  eigenfunction  expansion  on  a  finite  domain.  The  first  operators  are,  indeed  equivalent 
and  the  second  is  not  mathematically  equivalent  because,  it  cannot  derive  the  numerical  approximations.  The  Fractional 
Diffusion  equations,  for  numerical  methods  are: 
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0  <t<T,  0  <  x  <  1  <  a  <  2 

u(Oj  0  =  u(L,  t)  =  0, 
u(x,  0)  =  g(x), 

3.  THE  MATRIX  TRANSFORM  METHOD 

This  method  is  used  for  the  fractional  diffusion  equation,  to  propose  with  homogeneous  boundary  conditions.  It 
shows  that,  the  Matrix  transform  method  gives  the  correct  accurate  analytic  solution.  This  is  formed  by  orthonormal 


eigenfunctions,  from  which  the  Matrix  transform  method  and  analytic  solution  are  based  on.  In  this  subsection,  the  new 


technique  can  apply  to  the  fractional  diffusion  equation,  based  on  orthonormal  eigenfunctions.  The  matrix  is  also  a 
symmetric  positive  definite,  which  exist  here  as  a  nonsingular  matrix. 

Even  the  matrix  is  an  eigenvector,  symmetric  tridiagonal  matrix  and  eigenvalues  which  can  be  formed  in  an 
easiest  method.  Numerical  examples  with  a  finite  domain  are: 


0  <  t  <  T}  0  <  x  <n}  1  <  a  <2 


0)  =  2  ( nr  — 


u(Q}  t)  =  ufrr,  t)  =  0. 


n=  1 
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Figure  1 


In  this  figure  1,  the  analytic  solution  and  the  solutions  obtained  by  the  third  numerical  method  is  given.  From  this 
it  can  be  seen  that,  the  matrix  numerical  solution  is  in  good  agreement  with  the  analytic  solution. 
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4.  ANALYTIC  SOLUTION  FOR  MATRIX  TRANSFORM  METHOD 

The  Matrix  Transform  Method,  henceforth,  as  MTM  and  analytic  solution  are  based  on  orthonormal 
eigenfunctions,  but  only  the  maximum  errors  of  the  MTM  is  present  here.  The  backward  differentiation  formulas  are  used 
to  solve  the  ordinary  differential  equations,  by  using  transformed  fractional  equation. 

From  the  Table  1  and  Table  2,  it  shows  the  maximum  errors  of  the  MTM,  for  solving  the  fractional  partial 
differential  equation.  The  MTM  is  stable  and  convergent  for  solving  the  fractional  diffusion  equation,  based  on 
orthonormal  eigenfunctions.  The  analytical  solutions  of  the  fractional  partial  differential  equations  are  derived  by  LI  and 
L2  Matrix  Transform  Method,  which  is  difficult  to  evaluate. 


Table  1:  Errors  of  LI  MTM 


h=n/n 

MTM  LI 

7C/10 

2.217E-2 

71/20 

5.759E-3 

tt/40 

1.481E-3 

ti/80 

3.727E-4 

Table  2:  Errors  of  L2  MTM 


h=7r/n 

MTM  L2 

7C/10 

1.9963E-2 

71/20 

5.1683E-3 

tt/40 

1.3386E-3 

ti/80 

3.3519E-4 

Figure  2 

The  numerical  examples  of  fractional  partial  differential  equation  on  a  finite  domain  are: 


3u(x>  0 
dt 


=  Kt 


ad\x\ 


;u{x,  t), 


0  <t<T,  0<  x<n, 


1  <  a  <2 


u{Q}  t)  —  sin4x} 
u(0j  0  =  u(n}  t)  =  0. 


Impact  Factor  (JCC):  5.9876 


NAAS  Rating:  3.76 


Fractional  Partial  Differential  Equations  -  A  Study  by  Numerical  Methods 


21 


By  demonstrating  the  calculation  of  these  two  methods,  we  can  consider  the  following  Table  3,  for  the  maximum 
errors  of  the  standard  method,  to  solve  the  fractional  diffusion  equation  and  the  fractional  dispersion  equation. 

Table  3:  Errors  of  Standard  Method 


H=l/N 

A1/A2 

Standard 

Method 

1/50 

1.8144E- 

2 

2.8191E-3 

1/100 

9.4917E- 

3 

1.5093E-3 

1/200 

4.8584E- 

3 

7.7821E-4 

1/400 

2.4586E- 

3 

3.9459E-4 

The  impact  of  fractional  order  is  again  demonstrated,  which  is  the  example  of  the  RFADE,  with  a  different 
condition. 
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Figure  3 

Figure  3,  is  a  comparison  of  the  fractional  diffusion  equation  and  the  fractional  dispersion  equation,  with  a  finite 

domain. 

5.  CONCLUSIONS 

This  paper  concludes  that,  the  third  effective  numerical  method  to  solve  the  fractional  diffusion  equation  and  the 
fractional  dispersion  equation,  on  a  finite  domain,  with  orthonormal  Eigen  functions,  has  been  described.  Two  methods, 
which  were  defined  by  the  fractions,  were  discussed.  The  Numerical  method’s  results  have  demonstrated  the  effectiveness 
and  convergence  of  the  Matrix  Transform  methods.  The  methods  and  techniques,  which  were  discussed  in  this  paper,  can 
be  easily  applied  to  solve  the  fractional  partial  differential  equations,  even  in  other  kinds. 
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